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Transition rates in continuously driven steady states (relevant to sheared complex fluids) were 
derived in P, 0] by demanding that no information other than the microscopic laws of motion and the 
macroscopic observables of the system be used to describe it. This implies that the {noneqmUbrium) 
reservoir, to which the system is weakly coupled, is fully characterised by its mean energy and mean 
flux. While we expect the resulting prescription for the rates in continuous- and discretised-time 
models to be equivalent, it is not trivial to see this from the expression for the rates derived in 
We demonstrate this equivalence for a model of activated processes solved previously for continuous 
time [3] , thus demonstrating consistency of the theory . 

The rates of stochastic transitions between microstates of a complex fluid under shear, subject to non-equilibrium 
noise from a reservoir of the sheared fluid, were found in [3, Q to respect certain detailed-balance-like rules. The rules 
were applied to a toy model, described in 1], in which the fluid's discrete state-space has a comb-shaped connectivity of 
permitted transitions, as defined in Fig.[T] This models a system whose energetically favoured states are incompatible 
with a large flux. Whichever state is currently occupied by the system (that might be composed of many interacting 
particles), can be represented by a filled-circle, that has the appearance of a single particle hopping in a real-space 
landscape. The stochastic transitions occur with a certain probability in one time step. 

As in [H (but now for discrete time), we define an equilibrium version of the model, and derive a driven-ensemble 
version. In the latter case, the non-equilibrium constraint is to demand that, after a long time r, the ensemble- 
averaged strain has a large positive value 70, defining a shear current 7 by 70 = jr. As shown in Fig. [1] the system 
can occupy states of two different types: high-energy states a that are connected to other a states, and low-energy 
traps p. It may seem that one is free to choose independently four parameters that define the model: the probabilities 
u, d, I, r of transitions upward, downward, leftward and rightward in Fig. [1] respectively. At equilibrium, detailed 
balance requires F'^ = r'"^ = w and it®'^ = d^'^ ex-p{—E) = p with the thermal energy scale defined by fcsT = 1, so that 
only two free parameters remain. Transition probabilities are therefore similarly restricted in the driven steady-state 
ensemble [ij , where a transition from microstate a to & has probability per discrete time-step 



Pal = pII exp [vjab + qb{v) - qa{v) - Q{v)^t] . 



(1) 



Here p'^ is the probability at equilibrium for the same transition, 7^^ the net flux gained by it, is a parameter 
conjugate to the flux, and qa{v) defined for discrete transitions by 
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a quantity measuring the system's propensity to exhibit flux in future from the initial state a. Q{v) is a property of 
the system's steady state and is related to the average current by dQ/dv = 7. Here, the equilibrium Green's function 
Ga{l, t), is the probability that the system will undergo a total strain 7 in time r, beginning from state a. For discrete 
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FIG. 1: The discrete comb-shaped state-space for a simple model system studied in jl| for continuous time, and here for discrete 
time. 
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time steps {At = 1) Eq. ^ yields for the four transitions probabilities: 

u = pe«°(''^-«''('')-'3M . (3) 

To evaluate these we must find qa{v) — and Q{v). 

The connectivity of the state-space (Fig. [1]) implies relationships between the various Green's functions (simplified 
in Q but given here in full). A system initially in state /3 hops to state a at time t with probability p{t) = p{l —pY~^ 
and is governed by the propagator Gaij, t) from there on. This translates to 

r r 

Gpil, r)^J2 G.(7, T - <) = ^ p (1 - G„(7, r-t). (4) 
t=i t=i 

In terms of the quantity mi(zy, t) = In X]^-oo ^iil^ ''') ^^'^ i defined for each of the states (i = a, /?), Eq. (j4|) yields 
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jm0(i/,r+l)-m;3(i/,r) _ ^ _|_ ^ 



In the limit r — > c» , rafi{v^ t + 1) — rafiiy, r) — > (5(i') and mi{h', t) — mj{v, t) '?i(j^) — ■ Hence, 



^<lc,(v)-qfj(v) _ 1 
P 



(5) 



(6) 



Since only the a states carry flux, the mean shear flux 7 — {r — I) fa = (r — I) u/(u + d) is determined solely by r and 
I and the occupancy fa of the a state. Using dQ/dv — 7, this gives a differential equation for Q : 

dQ 2uj sinhz^e'^ 

which can be solved for Q{v). The condition Q(0) = (coming from normalization of the propagator) sets the 
constant of integration, yielding 



= uj cosh{v) + 1 - - |(e^ + 1) + [uocosh{v) + ^(1 - e-^) - cjj + p2 (g) 

Eliminating in favour of 7, Eq. thus gives the following four simultaneously equations that uniquely specify the 
rates in the driven steady state; 

r-l = {l + d/u)^ (9) 
[du) = p^e^ {rl) (10) 
u = l-\/W(l-p)/w (11) 
u;{r + l + d-u) = VTl [2 uj + p {e^ - 1)] . (12) 

The quantities r, Z, d, u are dimcnsionless, normalized probabilities, in contrast to the dimensionful rates analysed 
in the continuous-time version of the model [l| . In the limit of short time-steps (where the unit of time is redefined 
to be At — > 0), all hopping probabilities per time-step vanish in proportion to At. We define the transition rates to 
be the corresponding probabilities scaled by 1/At (denoted by i?, L, D, U as in 1]) with units chosen such that 
= L°'i = 1 in the continuous-time version of the model. The rates are plotted in Fig. [2^ for parameter values 
consistent with Fig. (3) of [1]. As At 0, our rates converge to those calculated previously for continuous time. 
In addition, it is easy to show that the constraints Eqs.([9l fT2|) lead to those in fH exactly, in the limit At — ^ , thus 
confirming the self-consistency of the formalism in [1| . 

The constraint of normalization in the discrete-time case leads to significant quantitative differences in the model's 
behaviour for non-vanishing time steps in the limit of large flux. The probabilities are plotted again in Fig[2}3 on 
logarithmic axes (for positive velocities), where the parameters E — 6, uj = 0.001, pe^ = 0.1 (as in [l|) have been 
chosen to separate the relevant time-scales for clarity. Vertical lines on the figure separate three different regimes of 
imposed current 7 (as also seen in the continuous-time version of the model [Ij): (a) the near-equilibrium regime at 
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FIG. 2: (a) The transition rates - probabilities for the discrete-time model scaled by 1/At - plotted for parameter values 
E = 2, pe^/co — 0.5 for two values of uj — 0.4 At) (dashed lines) and uj = 0.04 = At (dotted lines) along with the rates 
of the continuous time version (solid lines), (b) Log- log graphs of u, d, I and r for discrete-time hopping in the comb-shaped 
state-space of Fig. [T] as a function of mean shear flux 7, with E = 6, lo = 0.001, pe^ = 0.1. The dotted lines separate three 
distinct regimes of behaviour. 

the lowest velocities where the rates of activation and relaxation continue to be governed by the equihbrium principle 
of detailed balance, (b) a regime where ^ 7 ^ pe^, in which r and I remain approximately constant with increasing 
7, while the required flux is achieved by variations of only the activation and relaxation probabilities u and d, and (c) a 
regime where all four probabilities vary — unlike the continuous-time result, the flux-inhibiting transition probabilities 
I and d strictly vanish as 7 — > while the favoured probabilities r and u rise towards unity. 

In summary, we have applied the rules that were derived for transition rates in [J] for a particular class of nonequilib- 
rium steady states to a simple model system with a comb-shaped state-space which it explores in discrete time-steps. 
Even for such a simple system, the application of the transition-rate rules is a non-trivial procedure. (See [2] for some 
recent simplifications.) Having carried out that procedure, we have been able to show that in the limit of small time 
steps the results tend to those derived for the continuous time version of the same model thereby demonstrating the 
self-consistency of the formalism. 

Acknowledgments This work was funded by the Royal Society and by EPSRC grant GR/T24593/01. 



[1] Evans RML, Detailed balance has a counterpart in non-equilibrium steady states, 2005 J. Phys. A: Math. Gen. 38, 293. 
[2] Evans RML, Rules for Transition Rates in Nonequilibrium Steady States, 2004 Phys. Rev. Lett. 92, 150601. 
[3] Baule A and Evans RML, Invariant quantities in shear flow, 2008 Phys. Rev. Lett. 101, 240601. 



